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ABSTRACT 

Unprecedented  observations  of  Hurricane  Isabel  (2003)  at  category  5  intensity  were  collected  from  12  to 
14  September.  This  study  presents  a  detailed  analysis  of  the  inner-core  structure,  atmospheric  boundary 
layer,  sea  surface  temperature,  and  outflow  layer  of  a  superintense  tropical  cyclone  using  high-resolution  in 
situ  flight-level,  NCAR  GPS  dropwindsonde,  Doppler  radar,  and  satellite  measurements.  The  analysis  of 
the  dropwindsonde  and  in  situ  data  includes  a  comprehensive  discussion  of  the  uncertainties  associated  with 
this  observational  dataset  and  provides  an  estimate  of  the  storm-relative  axisymmetric  inner-core  structure 
using  Barnes  objective  analysis.  An  assessment  of  gradient  and  thermal  wind  balance  in  the  inner  core  is 
also  presented.  The  axisymmetric  data  composites  presented  in  this  study  suggest  that  Isabel  built  a  res¬ 
ervoir  of  high  moist  entropy  air  by  sea-to-air  latent  heat  flux  inside  the  low-level  eye  that  was  utilized  as  an 
additional  energy  source  to  nearly  maintain  its  extreme  intensity  even  after  crossing  the  cool  wake  of 
Hurricane  Fabian.  It  is  argued  here  that  the  combined  mean  and  asymmetric  eddy  flux  of  high  moist  entropy 
air  from  the  low-level  eye  into  the  eyewall  represents  an  additional  power  source  or  “turbo  boost”  to  the 
hurricane  heat  engine.  Recent  estimates  of  the  ratio  of  sea-to-air  enthalpy  and  momentum  exchange  at  high 
wind  speeds  are  used  to  suggest  that  Isabel  utilized  this  extra  power  to  exceed  the  previously  assumed 
intensity  upper  bound  for  the  given  environmental  conditions  on  all  three  days.  This  discrepancy  between 
a  priori  potential  intensity  theory  and  observations  may  be  as  high  as  35  m  s_1  on  13  September. 


1.  Introduction 

Category  5  and  supertyphoon-class  tropical  cyclones1 
(TCs)  are  some  of  the  most  awe-inspiring  natural  phe¬ 
nomena  on  the  planet.  Their  savage  beauty  and  tre¬ 
mendous  power  have  had  profound  impacts  on  history 
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(Emanuel  2005),  and  these  relatively  rare  phenomena 
have  the  potential  damage  value  of  500  times  that  of  a 
category  1  storm  (Pielke  and  Landsea  1998).  Meteoro¬ 
logical  observations  from  storms  that  achieve  this  infa¬ 
mous  status  provide  a  glimpse  of  the  hurricane  heat 
engine  operating  at  peak  efficiency,  yielding  new  in¬ 
sights  into  the  dynamics  and  thermodynamics  of  TCs. 

One  goal  toward  reaching  a  complete  understanding 
of  hurricane  intensity  is  developing  an  accurate  theory 
that  predicts  a  reasonable  upper  limit,  or  the  potential 
intensity  (PI),  of  a  TC  for  a  given  set  of  environmental 
conditions.  PI  theory  not  only  provides  a  prediction  of 


1  Category  5  refers  to  the  highest  level  on  the  Saffir-Simpson 
scale  in  the  Atlantic  basin,  while  supertyphoon  is  a  term  used  for 
intense  typhoons  in  the  western  Pacific  basin.  They  are  defined  by 
maximum  sustained  1-min  surface  winds  greater  than  135  and  130 
kt,  respectively. 
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the  capability  for  a  storm  to  achieve  category  5,  but  also 
a  simplified  framework  in  which  to  study  the  processes 
that  compose  the  hurricane  engine.  Understanding  in¬ 
tensity  change  then  becomes  a  matter  of  accurately  de¬ 
termining  the  limiting  factors  that  prevent  a  TC  from 
realizing  and/or  maintaining  its  peak  intensity  and  effi¬ 
ciency. 

One  approach  to  PI  is  to  use  statistical  predictors, 
such  as  the  work  of  DeMaria  and  Kaplan  (1994).  As  a 
testament  to  the  difficulty  in  forecasting  TC  intensity, 
statistical  models  currently  show  more  skill  than  opera¬ 
tional  dynamical  models  (DeMaria  et  al.  2005).  Statis¬ 
tical  PI  provides  an  important  empirical  upper  limit  for 
intensity  and  verification  of  the  importance  of  environ¬ 
mental  parameters  that  are  favorable  for  development, 
but  it  has  little  to  say  about  the  processes  that  control 
TC  intensity.  In  an  effort  to  understand  these  funda¬ 
mental  mechanisms,  several  other  approaches  have 
been  taken  over  the  last  few  decades  to  predict  the  PI 
based  on  simplifications  of  the  energetics  and  dynamics 
of  the  atmosphere-ocean  system.  These  theories  should 
be  valid  for  a  range  of  environmental  conditions,  and 
significant  violations  of  the  predicted  PI  beyond  obser¬ 
vational  error  are  therefore  not  statistical  anomalies, 
but  indications  that  the  dynamical  basis  of  the  theory  is 
either  flawed  or  incomplete. 

Predictions  for  the  minimum  surface  pressure  using  a 
parcel  approach  were  first  developed  by  Miller  (1958). 
The  surface  pressure  was  obtained  through  two  hydro¬ 
static  integrations:  by  first  lifting  an  air  parcel  moist 
adiabatically  in  the  eyewall,  and  subsequently  sinking 
dry  adiabatically  (with  mixing  from  the  eyewall)  in  the 
eye.  Holland  (1997)  modified  this  theory  by  incorpo¬ 
rating  the  pressure  dependence  of  moist  entropy  and 
reducing  the  surface  pressure  in  the  eyewall  integration 
iteratively  until  the  solution  converged.  Neither  of 
these  two  approaches  explicitly  considered  heat  and 
momentum  transfer  between  the  atmosphere  and 
ocean.  In  contrast,  heat  and  momentum  transfer  plays  a 
primary  role  in  theories  proposed  by  Kleinschmidt 
(1951),  Malkus  and  Riehl  (1960),  and  Emanuel  and  col¬ 
leagues  (Rotunno  and  Emanuel  1987;  Emanuel  1986, 
1988,  1991,  1995,  1997;  Bister  and  Emanuel  1998; 
Emanuel  et  al.  2004;  this  theoretical  work  is  hereinafter 
collectively  referred  to  as  EPI).  We  have  focused  this 
paper  on  comparisons  with  EPI  theory  since  air-sea 
heat  and  momentum  fluxes  are  not  explicitly  included 
in  Miller’s  (1958)  or  Holland’s  (1997)  PI  theories,  nor 
do  these  theories  offer  predictions  for  maximum  sus¬ 
tained  tangential  wind  (Umax).  Camp  (1999)  and  Camp 
and  Montgomery  (2001)  provide  a  more  comprehen¬ 
sive  review  of  PI  theory  for  interested  readers. 


Previous  studies  have  tested  various  PI  theories  using 
numerical  models  (Rotunno  and  Emanuel  1987;  Pers- 
ing  and  Montgomery  2003,  hereinafter  PM03),  east  Pa¬ 
cific  Ocean  environmental  soundings  (Hobgood  2003), 
and  best-track  and  satellite  datasets  (Tonkin  et  al. 
2000).  PM03  performed  a  rigorous  test  of  EPI  theory 
with  the  axisymmetric  numerical  model  developed  by 
Rotunno  and  Emanuel  (1987)  and  found  that  the  the¬ 
oretical  PI  was  exceeded  when  the  resolution  of  the 
model  was  able  to  resolve  the  hurricane  eye  and  its 
interaction  with  the  eyewall.  Their  work  suggested  that 
the  low-level  eye,  traditionally  passive  in  PI  theories 
(e.g.,  Emanuel  1995),  plays  an  important  role  in  the 
energetics  of  the  TC  engine  by  providing  a  reservoir  of 
high  moist  entropy  air  that  augments  the  energy  avail¬ 
able  between  the  outer  core  and  the  base  of  the  eye- 
wall.  PM03  coined  the  term  “superintensity”  to  refer  to 
storms  that  exceeded  their  PI  by  utilizing  this  mecha¬ 
nism. 

Comparable  tests  of  PI  with  observational  data  are 
limited.  A  climatological  study  of  hurricane  intensity  by 
Tonkin  et  al.  (2000)  showed  that  both  the  Emanuel  and 
Holland  PI  theories  demonstrated  reasonable  predic¬ 
tive  capability  for  all  ocean  basins  and  seasons,  with  a 
general  tendency  for  overestimation  of  the  intensity. 
Given  the  wide  range  of  adverse  environmental  condi¬ 
tions  a  tropical  cyclone  might  encounter,  this  result  is 
perhaps  not  too  surprising.  However,  an  interesting 
subset  of  storms  was  found  that  exceeded  their  PI. 
Their  study  showed  a  broad  range  of  SSTs  for  which 
both  PI  models  underestimated  the  TC  intensity. 
Whether  this  is  due  to  the  use  of  climatological  sea 
surface  temperature  (SST),  satellite-based  intensity  es¬ 
timates,  or  a  flawed  theory  is  not  known.  More  detailed 
observations  are  therefore  needed  to  elucidate  the 
structure,  intensity,  and  environmental  parameters  of 
the  anomalous  cases. 

A  critical  limitation  to  testing  the  limits  of  PI  theory 
is  the  lack  of  detailed  observations  of  the  inner-core 
kinematic  and  thermodynamic  structure  and  environ¬ 
mental  conditions  of  category  5  TCs.  Even  with  opera¬ 
tional  aircraft  reconnaissance,  it  is  difficult  to  establish 
high-resolution  two-  or  three-dimensional  structures  of 
a  hurricane  (Hawkins  and  Rubsam  1968;  Hawkins  and 
Imbembo  1976;  Frank  1984;  Jorgensen  1984;  LeeJoice 
2000). 2  Intensive  observing  periods  (IOPs)  conducted 
as  part  of  a  dedicated  field  project  are  an  effective,  and 
often  only,  way  to  obtain  the  quantity  and  quality  of 


2  Airborne  Doppler  radars  can  now  provide  operational  esti¬ 
mates  of  three-dimensional  kinematic  and  precipitation  fields,  but 
high-resolution  multidimensional  spatial  fields  of  thermodynamic 
data  are  still  difficult  to  obtain. 
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observations  needed  to  adequately  test  meteorological 
theory.  Three  IOPs  into  Hurricane  Isabel  were  ob¬ 
tained  from  12  to  14  September  2003  while  it  was  at 
category  5  intensity.  These  IOPs  were  supported  by  the 
National  Oceanic  and  Atmospheric  Administration 
(NOAA)  and  Office  of  Naval  Research  (ONR)  as  part 
of  the  Coupled  Boundary  Layer  Air-Sea  Transfer 
(CBLAST)  field  campaign,  NOAA  Hurricane  Field 
Program,  and  NOAA/National  Environmental  Satellite, 
Data,  and  Information  Service  (NESDIS)  Ocean 
Winds  experiment. 

The  Isabel  dataset  provides  a  unique  opportunity  to 
use  detailed  in  situ  observations  of  the  inner-core  struc¬ 
ture,  atmospheric  boundary  layer,  sea  surface  tempera¬ 
ture,  and  outflow  layer  to  test  the  predictions  of  Eman¬ 
uel’s  PI  theory,  and  examine  the  structure  and  balance 
of  the  hurricane  engine  operating  at  near-peak  effi¬ 
ciency.  This  study  will  quantify  the  environmental  pa¬ 
rameters  used  in  the  EPI  theory  using  measurements 
obtained  during  the  CBLAST  field  campaign,  and  com¬ 
pare  the  resulting  maximum  intensity  predictions  with 
the  observed  intensity  of  a  category  5  TC.  These  data 
are  also  used  to  assess  the  validity  of  the  assumptions  of 
gradient  and  thermal  wind  balance  invoked  in  the  der¬ 
ivation  of  the  theory.  EPI  theory  predicts  an  upper 
bound  on  the  mean  tangential  wind  at  the  swirling 
boundary  layer  top  for  given  environmental  conditions. 
Montgomery  et  al.  (2006,  hereafter  M06)  suggested  that 
Hurricane  Isabel  surpassed  this  upper  bound  on  13  Sep¬ 
tember  2003,  providing  the  first  evidence  for  the  exis¬ 
tence  of  superintense  storms  in  nature  that  were  previ¬ 
ously  predicted  using  high-resolution  computer  simula¬ 
tions  (PM03).  This  paper  extends  this  analysis  over 
three  consecutive  days  from  12  to  14  September  2003 
and  suggests  that  Isabel  was  superintense  over  the  en¬ 
tire  period. 

The  in  situ  fields  of  pseudoequivalent  potential  tem¬ 
perature  (labeled  hereafter  de )  and  radial  wind  ana¬ 
lyzed  here  suggest  that  a  reservoir  of  high  entropy  air 
was  enhanced  through  sea-to-air  latent  heat  flux  in  the 
low-level  eye  between  12  and  13  December,  and  that 
significant  penetration  of  near-surface  air  from  the  in¬ 
flow  was  enhanced  thermodynamically,  acquiring  the 
characteristics  of  the  high  entropy  air  in  the  eye.  This 
low-level  inflow  replenishes  the  air  transported/mixed 
out  of  the  eye  and  provides  additional  power  to  the 
hurricane  by  injection  into  the  eyewall  cloud.  The  eye 
dynamics  can  therefore  be  figuratively  compared  to  a 
second  cycle  or  “afterburner”  for  the  Carnot  engine,  in 
which  thermodynamic  energy  drawn  from  the  underly¬ 
ing  ocean  within  the  eye  augments  the  energy  obtained 
from  the  ocean  underneath  and  outside  the  eyewall, 
where  current  theory  assumes  all  of  the  energy  uptake 


occurs  (e.g.,  Emanuel  1995,  1997).  It  is  suggested  that 
the  low-level  eye  was  utilized  as  an  additional  energy 
source  to  nearly  maintain  its  extreme  intensity  through 
14  September,  even  after  crossing  the  cool  wake  of 
Hurricane  Fabian. 

Section  2  details  the  dataset  and  analysis  methodol¬ 
ogy  used  in  this  study.  The  evolution  of  Isabel’s  axi- 
symmetric  structure  from  12  to  14  September  is  pre¬ 
sented  in  section  3,  including  what  we  believe  is  the 
highest-resolution  potential  vorticity  fields  of  the  hur¬ 
ricane  inner  core  derived  from  observational  data  to 
date.  An  analysis  of  gradient  and  thermal  wind  balance 
is  given  in  section  4.  Observed  environmental  condi¬ 
tions  and  a  summary  of  EPI  predictions  are  presented 
in  section  5.  A  summary  and  conclusions  are  given  in 
section  6. 

2.  Dataset  and  analysis  methodology 

a.  Hurricane  Isabel  (2003) 

Hurricane  Isabel  became  a  tropical  depression  from 
an  African  easterly  wave  at  0000  UTC  6  September 
2003  and  was  quickly  named  a  tropical  storm  6  h  later 
(Lawrence  et  al.  2005).  Isabel  continued  intensification 
until  11  September,  when  Dvorak  satellite  estimates 
indicated  it  reached  category  5  intensity,  with  an  esti¬ 
mated  minimum  central  pressure  of  915  hPa  and  maxi¬ 
mum  sustained  surface  winds  of  145  kt.  Isabel  main¬ 
tained  maximum  sustained  surface  wind  speeds  above 
130  kt  and  central  pressures  below  940  hPa  in  relatively 
favorable  environmental  conditions  until  15  Septem¬ 
ber,  at  which  point  environmental  vertical  wind  shear 
increased  and  the  storm  began  to  weaken.  Isabel  made 
landfall  in  North  Carolina  three  days  later  as  a  large 
category  2  hurricane.  The  National  Hurricane  Center 
(NHC)  best  track  and  intensity  are  shown  in  Fig.  1.  Two 
NOAA  WP-3Ds  (P3s)  conducted  three  IOPs  from  12  to 
14  September  as  part  of  the  CBLAST  and  NOAA/ 
NESDIS  Ocean  Winds  experiments.  NHC  best-track 
intensities  were  estimated  at  135-140  kt  (category  5) 
during  the  —1600-2300  UTC  time  period  when  obser¬ 
vations  were  collected  on  each  day,  and  an  estimated 
130  kt  at  the  end  of  the  final  IOP  on  the  14th.  The 
85-GHz  and  visible  satellite  imagery  from  the  three 
analysis  days  are  shown  in  Fig.  2.  The  nearly  circular 
core  of  deep  convection  of  Isabel,  represented  by  the 
cold  (red)  85-GHz  brightness  temperatures  and  asym¬ 
metric  structures  in  the  low-level  stratus  clouds  in  the 
eye,  are  evident  during  all  three  IOPs.  While  the  mag¬ 
nitude  of  the  asymmetries  is  generally  an  order  of  mag¬ 
nitude  smaller  than  the  primary  circulation  in  the  inner 
core  (Shapiro  and  Montgomery  1993;  Reasor  et  al. 
2000),  the  low-level  vortex  structures  evident  in  the 
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Fio.  1.  NOAA  Tropical  Prediction  Center/National  Hurricane 
Center  (a)  best  track  and  (b)  best-track  intensity  for  Hurricane 
Isabel.  Open  hurricane  symbol  indicates  tropical  storm  strength, 
filled  symbol  indicates  hurricane  strength,  and  “L”  indicates  ex- 
tratropical  transition.  The  three  intensive  observing  periods  on  12, 
13,  and  14  Sep  are  highlighted. 


data,  but  radar  analysis  was  also  performed  for  inde¬ 
pendent  verification  of  the  large-scale  wind  fields  (not 
shown)  and  small-scale  features  (Aberson  et  al.  2006). 
Pseudoequivalent  potential  temperature  was  calculated 
following  the  empirical  formulation  in  Bolton  (1980).3 

b.  In  situ  flight-level  data 

Flight-level  in  situ  data  used  in  this  study  were  kindly 
provided  by  NOAA’s  Flurricane  Research  Division. 
Data  were  available  at  1-s  resolution  for  NOAA  air¬ 
craft  and  10-s  resolution  for  the  Air  Force  C130s.  G-IV 
flight  level  data  were  not  used,  except  for  dropwind- 
sonde  release  information.  A  rudimentary  correction 
for  instrument-wetting  errors  (Zipser  et  al.  1981;  Eastin 
et  al.  2002)  was  applied  to  supersaturated  dewpoint 
temperature  measurements.  This  correction  assumes 
that  the  errors  for  each  sensor  are  equal  in  magnitude 
but  opposite  in  sign  and  was  shown  by  Eastin  et  al. 
(2002)  to  reduce  the  majority  of  significant  wetting  er¬ 
rors  but  not  remove  them  completely,  resulting  in  a 
mean  6e  error  of  2.7  K. 

c.  NCAR  GPS  dropwindsondes 

An  unprecedented  total  of  184  National  Center  for 
Atmospheric  Research  (NCAR)  global  positioning  sys¬ 
tem  (GPS)  dropwindsondes  (sondes)  were  released 
into  the  inner  and  outer  core4  of  Flurricane  Isabel  dur¬ 
ing  the  three  CBLAST  IOPs.  An  additional  38  sondes 
were  released  in  the  ambient  environment  by  the 
NOAA  G-IV  on  13  and  14  September.  This  instrument 
provides  pressure,  temperature,  relative  humidity 
(PTH),  and  horizontal  wind  speed  at  2-Hz  temporal 
resolution  along  a  Lagrangian  trajectory  falling  be¬ 
tween  12  and  15  ms-1  in  the  lower  troposphere.  This 
yields  a  vertical  resolution  of  approximately  5  m,  with 
typical  PTH  errors  less  than  1.0  hPa,  0.2°C,  and  5% 
respectively,  and  wind  errors  less  than  2.0  m  s_1  (Hock 
and  Franklin  1999).  Vertical  velocity  was  obtained  from 
the  dropwindsondes  by  removing  the  estimated  termi- 


stratus  clouds  within  the  eye  (Fig.  2,  right  column)  are 
argued  here  to  be  important  elements  that  help  to 
maintain  category  5  intensity  over  these  three  days. 

The  NOAA  G-IV  and  U.S.  Air  Force  (USAF)  C130 
aircraft  also  collected  observations  during  the  IOPs  for 
operational  surveillance  and  reconnaissance.  In  situ 
flight-level  and  dropwindsonde  data  were  collected  by 
all  four  aircraft,  with  additional  Doppler  radar  and  ra¬ 
diometer  data  obtained  by  the  P3s  only.  This  study  fo¬ 
cuses  on  the  in  situ,  dropwindsonde,  and  radiometer 


3  Strictly  speaking,  EPI  theory  is  derived  using  reversible  ther¬ 
modynamics.  Because  of  the  lack  of  liquid  water  measurements 
from  the  dropwindsonde  data,  a  pseudoequivalent  potential  tem¬ 
perature  is  used  throughout  this  study.  Comparisons  including 
liquid  water  measurements  from  the  Johnson-Williams  liquid  wa¬ 
ter  probe  at  the  700-hPa  flight  level  indicate  that  this  approxima¬ 
tion  leads  to  a  linear  discrepancy  between  the  reversible  and  pseu¬ 
doequivalent  temperatures  of  —0.14%  per  g  kg-1  of  total  water 
mixing  ratio.  This  indicates  pseudoequivalent  potential  tempera¬ 
ture  is  —11  K  larger  at  —20  g  kg~\ 

4  Inner-core  observations  are  defined  in  this  paper  as  <60-km 
radius  from  the  center,  outer  core  as  150-250-km  radius. 
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Fig.  2.  (left)  Satellite  appearance  of  Hurricane  Isabel  at  85  GHz  (courtesy  NRL/Monterey)  and  (right)  visible  (courtesy  CIRA/CSU) 
wavelengths  during  each  IOP.  The  85-GHz  imagery  is  from  (a)  TMI  at  2126  UTC  12  Sep,  (c)  SSM/I  at  2218  UTC  13  Sep,  and  (e)  TMI 
at  2110  UTC  14  Sep.  Visible  images  are  from  GOES  super-rapid-scan  operations  at  1745  UTC  on  each  day. 
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nal  fall  speed  of  the  sondes  as  a  function  of  pressure. 
This  technique  has  been  shown  to  be  relatively  robust 
for  mesoscale  vertical  motions  in  hurricanes,  with  an 
estimated  error  of  0. 5-1.0  m  s-1  (Franklin  et  al.  2003). 

All  dropwindsondes  were  quality  controlled  to  re¬ 
move  noise  and  other  instrument  errors  with  either 
NCAR  Aspen  or  NOAA/Flurricane  Research  Division 
(HRD)  Editsonde  software.  For  this  study,  all  NOAA- 
released  sondes  and  most  USAF  sondes  were  kindly 
processed  by  FIRD  using  Editsonde,  with  some  addi¬ 
tional  USAF  sondes  processed  by  the  authors  using 
Aspen.  Even  though  both  of  these  programs  are  based 
on  the  same  quality  control  algorithms,  tests  were  per¬ 
formed  to  determine  the  differences,  if  any,  between 
the  two  software  packages.  Though  individual  pro¬ 
cessed  results  occasionally  exhibited  some  differences, 
only  minor  discrepancies  were  found  in  the  resulting 
axisymmetric  composites  after  averaging  the  large 
number  of  data  points. 

The  relative  humidity  for  a  particular  set  of  drop¬ 
windsondes  released  from  the  NOAA  43  P3  aircraft 
was  problematic,  however,  reaching  only  —85%  RF1  in 
the  eye  wall,  with  values  as  low  as  65%.  Since  the  the¬ 
oretical  value  in  precipitation  and  thick  marine  stra- 
tocumulus  should  be  near  100%,  this  suggested  a  po¬ 
tential  dry  bias.  Flowever,  the  true  homogeneity  of  the 
eyewall  cloud  is  unknown,  and  arbitrary  saturation  of 
eyewall  profiles  was  not  justified.  A  scaled  RF1  correc¬ 
tion  of  — 12%-15%,  was  therefore  applied  to  the  af¬ 
fected  dropwindsondes  (see  Bell  2006  for  more  details). 
The  correction  was  stable,  suggesting  that  the  degree  of 
molecular  contamination  (Wang  2005)  was  consistent 
for  that  batch  of  sondes.  A  few  minor  RFI  adjustments 
were  made  on  some  sondes  released  by  NOAA  42  as 
part  of  routine  quality  control.  While  the  applied  cor¬ 
rection  seems  plausible,  the  degree  to  which  the  cor¬ 
rected  sonde  profiles  represent  the  true  eyewall  RFI  is 
still  somewhat  uncertain.  To  verify  that  the  analysis  was 
not  overly  sensitive  to  the  uncertainty  of  these  mea¬ 
surements,  two  additional  composites  were  constructed 
without  the  suspect  data  and  without  applying  any  RFI 
corrections.  These  led  to  changes  of  de  in  the  eyewall  of 
2-4  K,  comparable  to  the  error  for  uncorrected  instru¬ 
ment  wetting  (Eastin  et  al.  2002).  It  is  therefore  be¬ 
lieved  that  the  postprocessing  of  the  relative  humidity 
data  removed  all  major  errors,  but  potential  sensor  wet¬ 
ting  and/or  molecular  contamination  still  yield  an  esti¬ 
mated  6e  uncertainty  of  —3  K. 

d.  Diagnosed  TC  centers 

Cylindrical  coordinates  allow  for  the  representation 
of  the  wind  field  by  tangential  and  radial  components, 
and  structural  decomposition  into  azimuthal  harmonics 


(wavenumbers).  This  requires  accurate  center  esti¬ 
mates  for  a  meaningful  coordinate  transformation.  The 
center-finding  method  of  Willoughby  and  Chelmow 
(1982,  hereafter  WC82)  relies  on  both  pressure  and 
wind  information  from  high-resolution  flight-level  data, 
and  has  been  shown  to  be  accurate  to  —3  km.  The 
variable  time  interval  between  fixes,  variability  of  the 
center  with  height,  and  local  pressure  and  wind  minima 
associated  with  mesovortices  in  the  eye  makes  deter¬ 
mining  an  accurate  set  of  centers  for  an  entire  IOP 
challenging  and  may  increase  the  center  uncertainty. 
The  analyzed  set  of  centers  for  each  day  was  created  by 
a  linear  interpolation  between  selected  robust  center 
fixes  obtained  by  the  WC82  method  using  storm- 
relative  winds  (i.e.,  winds  with  the  storm  motion  vector 
removed)  at  2-km  height,  with  constraints  provided  by 
NHC  best-track  data  when  reliable  center  fixes  were 
not  available.  This  yielded  a  general  west-northwesterly 
storm  motion  of  —7  m  s-1  on  each  day. 

Errors  in  the  analysis  introduced  by  center  uncer¬ 
tainty  were  examined  by  a  Monte  Carlo  approach,  per¬ 
turbing  the  estimated  center  to  account  for  center  un¬ 
certainties  resulting  from  small-scale  oscillations,  inter¬ 
polation  error,  and  vortex  tilt.  This  yielded  minor 
differences  for  reasonable  center  displacements  of  up 
to  5  km.  Dropwindsondes  released  in  the  eye  were  ex¬ 
amined  individually  for  potential  errors  in  center  posi¬ 
tion,  because  of  the  increased  sensitivity  of  the  cylin¬ 
drical  coordinate  transform  at  small  radii.  Given  the 
accuracy  of  the  individual  center  fixes  and  relative  ro¬ 
bustness  of  axisymmetric  quantities  to  the  random  per¬ 
turbation,  there  appear  to  be  no  systematic  errors  in¬ 
troduced  by  an  estimated  <5-km  mean  center  uncer¬ 
tainty. 

e.  Barnes  objective  analysis 

After  quality  control  and  decomposition  into  storm- 
relative  cylindrical  coordinates,  the  resulting  data  dis¬ 
tribution  in  the  radial-azimuthal  and  radial-vertical 
planes  for  each  IOP  is  shown  in  Fig.  3.  By  making  the 
steady-state,  axisymmetric  assumption,  the  primarily 
Eulcrian  flight  level  and  Lagrangian  dropwindsonde 
measurements  can  both  be  treated  as  axisymmetric,  in¬ 
stantaneous  snapshots  in  the  radial-vertical  plane  of 
the  storm  circulation.  The  aircraft  or  dropwindsonde 
measurement  is  then  weighted  according  to  its  distance 
from  a  given  storm-relative  radius  and  altitude  (grid 
point).  The  weighted  observations  from  different  azi¬ 
muths  and  times  are  then  averaged  by  the  Barnes  ob¬ 
jective  analysis  (Barnes  1973;  Koch  et  al.  1983)  to  yield 
an  estimate  of  the  axisymmetric  structure.  Whereas  the 
Barnes  objective  analysis  has  been  used  previously  in 
hurricane  studies  (e.g.,  Velden  et  al.  1992;  Franklin  et 
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Fig.  3.  Dropwindsonde  locations  and  trajectories  and  aircraft  flight  tracks  relative  to  storm  center 
from  1600  to  2300  UTC  on  each  day.  Storm-relative  data  distribution  in  the  (left)  azimuthal  (r-0)  plane 
and  (right)  radial-height  (r-z)  plane,  showing  the  NOAA  P-3  (42  in  blue,  43  in  green),  USAF  C-130 
(in  black)  flight  tracks,  and  dropwindsonde  trajectories  (in  red).  The  dropwindsondes  in  the  left  column 
move  cyclonically  (counterclockwise). 
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al.  1993),  its  application  for  inferring  radius-height 
mean  hurricane  structure  is  believed  to  be  novel.  A 
total  of  —43  500,  30  700,  and  37  800  data  points  were 
available  on  each  day  from  dropwindsonde  and  flight- 
level  data.  The  sounding  data  comprised  63%,  48%, 
and  59%  of  the  data  distributions  for  each  IOP,  from 
67,  35,  and  58  inner-core  soundings,  respectively.  The 
Barnes  analysis  procedure  used  radial  and  vertical  grids 
of  2500  and  250  m,  respectively,  to  account  for  the  dif¬ 
fering  horizontal  and  vertical  data  resolution  and  spac¬ 
ing.  A  minimal  gamma  smoothing  parameter  of  0.3  and 
radial  and  vertical  weight  parameters  were  then  set  to 
produce  the  maximum  spectral  resolution  (10  and  1  km, 
respectively)  allowable  for  the  given  grid  spacing. 

All  variables  that  do  not  contain  derivatives  were 
calculated  from  the  observational  data  and  then  com¬ 
posited  using  the  Barnes  analysis.  Given  the  relatively 
dense  observations  and  high  spatial  resolution  of  the 
composite,  an  attempt  was  also  made  to  assess  kine¬ 
matic  and  thermodynamic  gradients  derived  from  the 
composite  in  the  inner  core.  The  vorticity  and  potential 
vorticity  (which  also  involves  the  potential  temperature 
gradient)  resulting  from  these  calculations  are  there¬ 
fore  presented  in  section  3,  with  the  caveat  that  some  of 
the  detailed  structure  and  exact  magnitude  may  be  grid 
dependent. 

Uncertainties  arising  from  undersampling,  the  loca¬ 
tion  of  the  circulation  center,  and  dropwindsonde  ter¬ 
minal  fall  speed  may  produce  unbalanced  divergence 
and  vertical  velocity  fields  at  any  particular  point,  re¬ 
sulting  in  a  nonzero  residual  in  the  mass  continuity 
equation.  Unfortunately,  the  resulting  calculation  was 
too  sensitive  to  quantitatively  determine  this  residual. 
The  finescale  composite-derived  gradients  of  pressure 
and  Qe  in  the  balance  diagnosis  (section  4)  were  suspect, 
and  derivatives  were  calculated  from  radially  and  ver¬ 
tically  binned  averages  instead  of  the  composite  data. 
This  technique  is  believed  to  more  accurately  capture 
the  magnitude  of  the  gradients  at  a  specific  vertical 
level,  with  the  cost  of  decreased  radial  spatial  resolu¬ 
tion,  which  is  appropriate  for  the  assessment  of  balance 
at  the  larger  storm  scale. 

3.  Axisymmetric  structure  and  evolution  from 
12  to  14  September 

a.  Kinematic  structure  and  evolution 

This  analysis  of  Hurricane  Isabel  suggests  that  de¬ 
spite  the  relatively  steady-state  intensity,  the  storm 
structure  slowly  evolved  from  12  to  14  September.  Fig¬ 
ure  4  shows  the  radius-height  composite  storm-relative 
tangential  wind  (color),  radial  wind  (contour),  and  sec¬ 
ondary  circulation  (vector)  in  meters  per  second.  The 


origin  (0,  0)  denotes  the  storm  center  at  the  ocean  sur¬ 
face.  Isabel  crossed  the  cool  wake  of  Fabian  on  13  Sep¬ 
tember,  resulting  in  a  1°-2°C  reduction  in  sea  surface 
temperature  (see  Fig.  9),  and  one  might  expect  the  in¬ 
tensity  to  have  decreased,  but  this  analysis  indicates 
otherwise.  The  mean  peak  tangential  wind  at  the  radius 
of  maximum  wind  (RMW)  remained  very  strong  on  all 
three  days,  with  the  core  region  of  maximum  tangential 
winds  decaying  from  —80  to  74  m  s-1,  rising  from 
— 500-m  to  1-km  altitude,  and  expanding  from  —25-  to 
50-km  radius.  The  low-level  radial  inflow  increased  in 
both  depth  and  intensity  from  12  to  13  September,  but 
then  weakened  again  on  14  September.  A  persistent 
region  of  —5-10  m  s”1  outflow  just  above  the  boundary 
layer  near  the  RMW  is  evident  on  all  three  days.  The 
derived  vertical  velocity  is  qualitatively  consistent  with 
the  radial  divergence,  showing  weak  vertical  motions 
inside  the  eye  and  a  maximum  updraft  nearly  collo¬ 
cated  with  the  RMW  on  each  day.  Lowest-level  (rep¬ 
resentative  of  —100  m)  radial  inflow  of  20  m  s-1  located 
at  25-km  radius  from  the  center  on  13  September  sug¬ 
gests  significant  penetration  of  air  from  the  eyewall  into 
the  eye.  This  measurement  of  strong  inflow  inside  the 
eye  appears  to  be  robust  (Bell  2006),  but  because  of 
limited  sampling  it  may  not  be  a  quantitatively  accurate 
depiction  of  the  axisymmetric  inflow  at  these  radii. 

There  is  the  possibility  that  the  expansion  of  the 
RMW  was  associated  with  an  eyewall  replacement 
cycle  from  12  to  13  September  (Kozich  2006).  Radar 
reflectivity  imagery  [not  shown;  reflectivity  composites 
are  available  from  NOAA/ Atlantic  Oceanographic  and 
Meteorological  Laboratory  (AOML)/HRD  at  www. 
aoml.noaa.gov/hrd]  and  microwave  imagery  (Fig.  2a) 
suggests  that  an  outer  rainband  began  to  encircle  the 
primary  eyewall  late  on  the  12th.  However,  only  limited 
flight-level  data  and  microwave  imagery  between  IOPs 
is  available,  making  it  difficult  to  determine  the  details 
of  an  eyewall  replacement  and/or  expansion.  The  tran¬ 
sition  between  12  and  13  September  occurred  with  a 
— 10-hPa  rise  in  central  pressure,  and  a  brief  weakening 
but  rapid  recovery  of  the  maximum  tangential  wind 
speed.  A  more  detailed  analysis  of  this  transition  is 
beyond  the  scope  of  this  study. 

b.  Equivalent  potential  temperature  and  absolute 
angular  momentum 

Figure  5  shows  the  radius-height  composite  0e 
(color),  specific  absolute  vertical  angular  momentum 
(contour),  and  transverse  secondary  circulation  (vec¬ 
tor).  The  specific  absolute  vertical  angular  momentum 
(absolute  circulation)  is  M  =  rv  +  Vifr2,  where  r  is  the 
radius  from  the  vortex  center,  v  is  the  storm-relative 
tangential  velocity,  and  f  is  the  Coriolis  parameter. 
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Fig.  4.  Radius-height  azimuthal  mean  storm-relative  tangential  wind  (color),  radial  wind  (contour),  and  the  secondary  circulation 
(vector)  in  m  s_1  derived  from  GPS  dropwindsonde  and  flight-level  data  from  (a)  12,  (b)  13,  and  (c)  14  Sep  2003. 


Some  of  the  most  dramatic  changes  in  the  inner-core 
storm  structure  are  illustrated  here,  with  a  distinct  in¬ 
crease  in  the  low-level  9e  in  the  eye  from  12  to  13  Sep¬ 
tember,  followed  by  an  increase  in  eyewall  0e  on  14 
September.  The  radial  <),,  gradient  is  generally  negative 
throughout  all  three  days,  except  for  very  near  the  cen¬ 
ter  on  the  12th.  On  the  14th,  it  appears  as  if  the  9e  has 
been  "mixed  out,”  with  relatively  lower  values  found  in 


the  eye,  and  an  increase  at  the  eyewall.  These  figures 
suggest  that  there  were  significant  changes  in  the  mean 
moist  entropy  structure  over  these  three  days.  Radial 
profiles  of  mean  potential  temperature  and  vapor  mix¬ 
ing  ratio  at  the  lowest  composite  level  are  shown  in  Fig. 
6.  These  radial  profiles  suggest  that  the  increase  in 
mean  9e  after  the  12th  was  primarily  due  to  increased 
low-level  moisture  and  occurred  despite  a  rise  in  central 


2032 


MONTHLY  WEATHER  REVIEW 


Volume  136 


0  5  10  15 

(c)  14  September 


25  30  35 

Radius(km) 


50  55  *  60 

10  m  s'1  I 

K  Figure  is  stretched  4x  in 

350  353  356  359  362  365  368  371  374  the  vertical  for  clarity 


Fig.  5.  Radius-height  azimuthal  mean  storm-relative  de  (color;  K);  absolute  angular  momentum  (contour;  m2  s  1  X  106);  and 
transverse  secondary  circulation  (vector)  from  (a)  12,  (b)  13,  and  (c)  14  Sep  2003. 
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pressure  of  —10  hPa.  This  supports  the  hypothesis  of 
persistent  latent  heat  flux  from  the  underlying  ocean 
inside  the  low-level  eye  and  possibly  radial  moisture 
flux  due  to  mixing  from  the  eyewall.  Whereas  a  detailed 
calculation  of  the  residence  time  of  air  parcels  in  Isa¬ 
bel’s  eye  is  beyond  the  scope  of  this  study,  Cram  et  al. 
(2007)  showed  that  residence  times  in  the  eye  for  sig¬ 
nificant  0e  gain  in  their  numerical  simulation  of  Hurri¬ 
cane  Bonnie  were  commonly  on  the  order  of  40-60  min 
and  as  low  as  15  min.  This  suggests  that  latent  heat 


fluxes  in  the  eye  can  modify  inflowing  air  parcels  on 
even  relatively  short  time  scales. 

The  concurrent  increase  of  the  moist  entropy,  verti¬ 
cal  velocity,  and  radial  inflow  from  12  to  13  September 
suggests  a  positive  feedback  between  these  fields.  The 
addition  of  higher  Qe  to  the  eyewall  cloud  enhances 
(locally)  buoyant  updrafts  (Braun  2002;  Eastin  et  al. 
2005;  Smith  et  al.  2005),  which  is  believed  to  support  a 
stronger,  deeper  radial  inflow.  This  may  allow 
Lagrangian  parcels  to  more  easily  breach  the  high  in- 
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Fig.  6.  Azimuthal  mean  radial  profiles  of  potential  temperature 
and  vapor  mixing  ratio  at  the  lowest  composite  level  (0-250  m, 
representative  of  —100  m)  from  12  to  14  Sep  2003. 


ertial  stability  and  Rossby  elasticity  (Cram  et  al.  2007) 
at  the  eyewall  and  be  enhanced  by  the  moist  entropy 
reservoir  in  the  eye  before  returning  to  the  eyewall 
cloud. 

A  key  ingredient  to  understanding  hurricane  maxi¬ 
mum  intensity  is  the  amount  of  moist  entropy  a  parcel 
is  able  to  gain  in  the  TC  boundary  layer.  Malkus  and 
Riehl  (1960)  posited  that  the  entropy  gain  occurs  on  the 
inward  spiral  from  the  environment  to  the  eyewall. 
Conversely,  Emanuel  (1986)  argued  that  convective 
downdrafts  in  the  rainbands  prevent  significant  mate¬ 
rial  increase  of  moist  entropy  for  inflowing  parcels, 
such  that  most  of  the  gain  occurs  directly  under  the 
eyewall.  PM03  explored  a  third  mechanism,  in  which 
persistent  enthalpy  fluxes  in  the  low-level  eye  provide 
an  important  additional  source  of  moist  entropy  for  the 
hurricane  engine.  As  a  proxy  for  the  Lagrangian  change 
in  0e  for  a  boundary  layer  air  parcel,  the  observed  val¬ 
ues  were  averaged  over  the  lowest  kilometer  in  radial 
bins  corresponding  to  the  outer  core,  eyewall  (at  the 
RMW),  and  eye.  These  values  are  shown  in  Table  1. 
The  outer-core  values  of  mean  boundary  layer  9e  were 
similar  for  each  IOP  (—350-352  K).  This  indicates  that 
the  gain  in  0e  for  a  hypothetical  inflowing  air  parcel  in 
the  boundary  layer  from  the  outer  core  to  the  eyewall 
was  the  largest  on  the  12  September  (— 11  K)  and  small¬ 
est  on  13  September  (—4  K).  An  additional  —4  and  —14 
K,  respectively,  was  available  for  parcels  that  were  able 
to  access  the  eye,  however.  Given  that  the  tangential 
velocity  remained  nearly  steady  state,  implying  that  the 


Table  1.  Mean  S0e  between  different  regions  of  the  storm  on 
the  three  analysis  days.  Values  are  averaged  over  the  lowest  ki¬ 
lometer  and  are  relative  to  values  in  the  outer  core  (—351  K  at 
— 250-km  radius).  Eyewall  values  are  from  10-km-wide  radial  bins 
at  25, 45,  and  50  km,  respectively,  and  eye  values  are  from  variable 
bins  inside  these  radii.  These  values  can  be  viewed  as  a  proxy  for 
the  Lagrangian  0e  increase  of  a  hypothetical  inflowing  air  parcel 
that  either  ascends  directly  in  the  eyewall  or  is  able  to  access  the 
eye  before  being  returned  to  the  eyewall. 


Mean  S0e 
(outer  core  — 
eyewall) 

Mean  S0C 
(eyewall  — 
eye) 

Mean  SSe 
(outer  core  - 
eye) 

12  Sep 

11  K 

4  K 

15  K 

13  Sep 

4  K 

14  K 

18  K 

14  Sep 

6  K 

7  K 

13  K 

frictional  dissipation  was  similar  on  each  day,  the  re¬ 
quired  moist  entropy  to  maintain  the  intensity  on  the 
13th  thus  most  likely  originated  in  the  eye.  It  is  also 
interesting  to  note  that  the  inner-core  SST  estimates 
were  the  lowest  on  this  day  (27.5°C  versus  28.5°-29°C; 
see  section  5a),  suggesting  that  the  reduced  latent  heat 
fluxes  underneath  the  eyewall  were  compensated  for  by 
radial  fluxes  from  the  eye.  Isabel  appears  to  have 
settled  into  an  intermediate  structure  on  the  14th,  with 
an  — 6-K  increase  from  the  outer  core  to  the  eyewall  and 
an  additional  —7  K  into  the  eye.  This  analysis  suggests 
that  inflowing  parcels  could  potentially  gain  on  average 
13-18  K  of  6e  on  all  three  days  via  different  thermody¬ 
namic  pathways. 


c.  Potential  vorticity  and  absolute  vertical  vorticity 

The  axisymmetric  dry  Rossby-Ertel  potential  vortic¬ 
ity  (PV;  Hoskins  et  al.  1985)  was  calculated  from  the 
composite  gradients  of  axisymmetric  tangential  velocity 
and  potential  temperature  in  cylindrical  coordinates  ac¬ 
cording  to  the  definition 


PV  =  a(£abs)  ■  V0  =  a 


d(rv)  100  dv  001 

r'dr  +  f  dz  dz  dr\  ’ 


(1) 

where  a  denotes  the  specific  volume,  £abs  denotes  the 
absolute  vorticity  vector, /denotes  the  Coriolis  param¬ 
eter,  v  denotes  the  tangential  velocity,  and  0  denotes 
the  potential  temperature.  There  is  only  a  slight  de¬ 
crease  in  the  vorticity  and  PV  of  the  vortex  each  day,  as 
shown  in  Fig.  7,  due  to  the  slowly  expanding  wind  field 
and  weakening  of  the  radial  gradient  of  tangential 
wind.  The  composite  vertical  vorticity  is  generally  posi¬ 
tive  everywhere  in  the  domain,  with  the  exception  of  a 
small  area  near  10-km  radius  on  14  September.  This  is 
believed  to  be  an  artifact  of  the  analysis  resulting  from 
the  lack  of  data,  large  eye,  and  weak  radial  gradient  of 
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Potential  Vorticity  (Color,  PVU  (10“6  K  kg'1  m2  s'1)) 
and  Absolute  Vertical  Vorticity  (1  O'3  s_1 ) 

4  (a)1^September 


10  15  20  25  30  35  40  45  50  55 


PVU 


0  10  20  30  40  50  60  70  80 

Fig.  7.  Radius-height  azimuthal  mean  storm-relative  Rossby-Ertel  potential  vorticity  [colors;  potential  vorticity  units  (PVU  =  10“ 
K  kg-1  m2  s-1)]  and  absolute  vertical  vorticity  (contours;  1CT3  s_1)  from  (a)  12,  (b)  13,  and  (c)  14  Sep  2003. 


tangential  wind  in  this  region.  For  consistency  with  the 
other  plots,  these  data  points  were  left  in  the  final 
analysis  despite  the  suspect  values  since  the  resulting 
remaining  potential  vorticity  structure  was  consistent 
with  the  other  days. 

The  "hollow  tower”  or  ring  structure  of  PV  at  ~l-km 
altitude  is  clearly  evident  on  all  three  days,  suggesting 


that  the  vortex  meets  the  necessary  condition  for  mixed 
barotropic/baroclinic  instability  (Montgomery  and  Sha¬ 
piro  1995;  Schubert  et  al.  1999).  A  combined  barotro¬ 
pic/baroclinic  instability  and  nonlinear  evolution  is  be¬ 
lieved  to  be  responsible  for  the  mesovortices  evident  in 
visual  satellite  (Fig.  2,  right  column)  and  photographic 
images  (not  shown)  of  coherent  vortex  structures  in  the 
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eye,  polygonal  eyewall  structures  observed  by  radar 
(M06,  their  Fig.  lc),  and  a  smaller-scale  eyewall  “miso- 
cyclone”  (Aberson  et  al.  2006).  A  detailed  stability 
analysis  that  includes  the  presence  of  moisture  is  a  topic 
of  some  interest  and  will  be  reported  in  due  course. 
These  mesovortices  emerge  as  the  vortex  sheet  on  the 
inside  of  the  eyewall  breaks  down.  This  breakdown 
process  enhances  radial  mixing  and  penetration  of  low- 
level  inflowing  air  into  the  eye  region  (Rotunno  1984; 
Emanuel  1997;  Kossin  and  Schubert  2001;  Montgomery 
et  al.  2002) 

There  are  two  regions  of  high  PV  evident  in  Fig.  7:  a 
lower,  outer  maximum  associated  with  the  strong  radial 
shear  at  the  inner  edge  of  the  eyewall,  and  an  upper, 
inner  maximum  associated  with  the  subsidence  inver¬ 
sion  in  the  eye.  The  eye  PV  maximum  at  ~3-km  alti¬ 
tude  corresponds  well  to  a  transition  region  with  a 
sharp  decrease  in  relative  humidity  values  (not  shown), 
and  weak  subsidence/near-zero  vertical  velocity,  remi¬ 
niscent  of  the  tropopause.  Inspection  of  individual 
dropsonde  profiles  indicates  that  the  height  of  the  in¬ 
version  was  not  homogeneous,  although  the  dataset 
suggests  it  was  between  2  and  4  km  on  average.  While 
some  of  the  finescale  detail  and  exact  magnitude  of  the 
PV  fields  is  most  likely  dependent  on  the  particular 
data  sampling  and  compositing  technique,  the  consis¬ 
tent  structure  inferred  on  all  three  days  provides  con¬ 
fidence  that  the  gross  features  of  this  important  quan¬ 
tity  are  effectively  captured  by  this  analysis.  This  struc¬ 
ture  is  qualitatively  consistent  with  full  physics 
numerical  simulations  (Chen  and  Yau  2001;  Wang  and 
Zhang  2003).  The  latter  study  shows  a  similar  feature  in 
the  low-level  eye  of  Hurricane  Andrew  (1992)  at  —2.5- 
km  altitude,  and  the  bowl-shaped  maximum  just  inside 
the  RMW  typically  associated  with  very  intense  storms. 
The  variability  of  the  eye  thermal  inversion  and  the 
dynamical  role  of  the  PV  maximum  associated  with  it 
are  interesting  topics  for  future  research. 

The  PV  presented  here  is  the  traditional  dry  Rossby- 
Ertel  PV,  calculated  by  the  inner  product  of  the  abso¬ 
lute  vorticity  vector  and  gradient  of  potential  tempera¬ 
ture.  It  is  important  to  note,  however,  that  part  of  the 
simplicity  of  the  EPI  framework  arises  through  an  as¬ 
sumption  of  zero  moist  saturated  PV,  in  which  the  po¬ 
tential  temperature  is  replaced  by  saturated  equivalent 
potential  temperature  (0*)  in  Eq.  (1).  In  a  zero  0J-PV 
vortex,  lines  of  constant  absolute  angular  momentum 
and  saturated  moist  entropy  are  congruent  above  the 
boundary  layer,  allowing  for  a  PV  inversion  through 
knowledge  of  the  boundary  conditions  alone.  In  the 
EPI  framework  the  complete  balanced  wind  and  ther¬ 
mal  fields  can  be  obtained  through  knowledge  of  the 


radial  structure  of  9e  in  the  boundary  and  outflow  layers 
and  the  assumption  of  saturation  throughout  the  hur¬ 
ricane  eyewall  in  the  middle  troposphere.  The  above 
analysis  suggests  that  while  the  M  and  0e  lines  are  in¬ 
deed  nearly  parallel  throughout  much  of  the  inner  core, 
there  are  significant  deviations  from  zero  0e-PV,  even 
in  the  nearly  saturated  eyewall,  on  all  three  days.  These 
deviations  appear  most  pronounced  in  the  lower  tropo¬ 
sphere.  This  is  consistent  with  radial  mixing  of  6e  and 
the  addition  of  heat  energy  from  the  eye  into  the  eye- 
wall. 


4.  Assessment  of  axisymmetric  balance 

To  derive  EPI  theory,  several  assumptions  about  the 
axisymmetric  balance  of  the  tropical  cyclone  need  to  be 
made,  including  hydrostatic,  gradient  wind,  and  ther¬ 
mal  wind  balance.  Before  attempting  to  make  predic¬ 
tions  using  this  theory,  gradient  wind  balance  is  revis¬ 
ited  (Willoughby  1990,  1991;  Gray  1991)  and  thermal 
wind  balance,  as  derived  by  Kleinschmidt  (1951)  and 
Emanuel  (1986),  is  examined  here  with  this  high- 
resolution  dataset. 


a.  Gradient  wind  balance 

The  azimuthally  averaged  radial  momentum  equa¬ 
tion  with  Boussinesq  approximation  in  cylindrical  co¬ 
ordinates  is  given  by 


dp 


-2 

V 


Du 


-  —  = - /F  +  —  +  F, 

pdr  \  r  )  Dt 

F  =  —(v'2)/r  +  d{ru'2)/rdr  +  d{u'w')/dz ,  and 
Dii/Dt  =  du/dt  +  udu/dr  +  wdii/dz, 


(2) 

(3) 

(4) 


where  F  represents  mean  eddy  flux  divergence;  Du/Dt 
is  the  material  derivative  of  the  average  radial  wind;  u, 
v,  and  w  are  the  cylindrical  velocity  components;  t  is 
time; /the  Coriolis  parameter  (assumed  constant);  p  the 
total  pressure,  and  p  the  total  density.  Here,  the  overbar 
represents  the  azimuthal  average  and  primes  represent 
perturbations  from  the  average.  If  the  mean  and  per¬ 
turbation  wind  field  is  known,  Eq.  (2)  can  be  integrated 
radially  to  obtain  the  estimated  contributions  from  each 
term  to  the  overall  mean  pressure  deficit.  Here  F  can¬ 
not  be  determined  with  the  current  observations.  How¬ 
ever,  it  can  be  estimated  by  a  residual  given  the  ob¬ 
served  pressure  gradient  and  mean  flow  quantities. 

The  results  of  such  an  integration  from  composite 
wind  data  on  13  September  are  shown  in  Fig.  8.  The 
mean  wind  is  assumed  steady  during  the  composite  pe¬ 
riod.  An  inward  integration  of  Eq.  (2)  using  the  com- 
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(a)  Calculated  vs.  observed  pressure  deficit  at  ~100  m 
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Fio.  8.  Calculated  pressure  deficit  obtained  by  integrating  the  radial  pressure  equation  [Eq.  (4.1)]  inward  from  55-km  radius  with 
(dark  gray  dashed  curve;  squares)  and  without  (black  solid  curve;  circles)  advection  terms  vs  the  pressure  deficit  observed  by  drop- 
windsonde  data  (light  gray  dotted  curve;  diamonds)  at  (a)  100-m  and  (b)  2-krn  altitude.  Third-degree  polynomial  fits  (light  gray  thin 
curve)  of  the  observed  pressure  gradient  with  r 2  values  of  (a)  0.98  and  (b)  0.99  are  also  shown. 


posite  mean  winds,  from  55-km  radius  with  zero  as  the 
outer  boundary  condition,  at  (a)  —100  m  (0-250-m 
composite  level)  and  (b)  2-km  altitude  is  compared  to 
the  observed  pressure  gradient.  The  integration  was 
performed  with  only  the  cydostrophic  and  Coriolis 
terms  (black  solid  curve)  and  also  including  the  advec¬ 
tion  terms  ( uduldr  +  wduldz,  gray  dashed  curve).  The 
observed  pressure  gradient  was  calculated  in  two  ways: 
by  averaging  dropwindsonde  data  in  10-km  radial  and 
50-m  vertical  bins  (dotted  curve),  and  from  a  ID  Barnes 
analysis  of  flight-level  data  at  2  km.  The  derived  pres¬ 
sure  field  from  the  2D  Barnes  analysis  was  sensitive  to 
data  gaps  because  of  the  strong  vertical  pressure  gradi¬ 
ent.  A  ID  Barnes  analysis  of  dropsonde  data  yields 
similar  results  to  the  binned  averages,  but  underesti¬ 
mated  the  steep  gradient  near  the  inner  edge  of  the 
eyewall  where  data  were  limited. 

Figure  8a  suggests  the  importance  of  the  advection 
terms  and  the  relatively  significant  contribution  from 
the  residual  terms  in  the  boundary  layer.  The  trans¬ 
verse  advection  and  eddy  flux  divergence  act  together 
to  oppose  the  cydostrophic  and  Coriolis  terms  and  re¬ 
duce  the  radial  pressure  gradient.  These  terms  increase 
the  central  pressure  by  —6  and  —8  hPa,  respectively. 
While  some  of  the  residual  may  be  due  to  uncertainties 
in  the  estimated  central  pressure,  the  combined  eddy 
terms  are  roughly  of  the  same  order  and  sign  as  the 
advection  terms  in  the  boundary  layer.  At  2-km  alti¬ 


tude,  the  calculated  and  observed  pressure  deficits 
agree  well,  consistent  with  a  decrease  in  radial  accel¬ 
erations  and  eddy  momentum  fluxes.  Comparisons  be¬ 
tween  the  dropsonde  and  flight-level-derived  pressure 
gradients  suggest  that  both  calculations  capture  the 
general  decreasing  trend  well,  though  the  binned  aver¬ 
aging  may  slightly  underestimate  the  gradient.  Similar 
results  from  12  and  14  September  (not  shown)  suggest 
that  the  vortex  is  roughly  in  gradient  balance  above  the 
boundary  layer  in  agreement  with  prior  analyses  of 
other  intense  storms  (e.g.,  Willoughby  1990). 

Alternately,  one  can  calculate  the  gradient  wind  from 
the  observed  pressure  field.  The  last  two  terms  in  (2) 
are  identically  zero  for  a  steady-state  vortex  in  pure 
gradient  wind  balance,  yielding  a  quadratic  equation 
that  can  be  solved  for  the  gradient  wind  (Vg)  for  a  given 
radial  pressure  distribution.  To  compute  the  gradient 
wind,  the  binned  pressure  data  was  fit  to  a  third-degree 
polynomial  (Figs.  8a,b,  thin  solid  line)  and  an  analytical 
derivative  was  obtained.  This  derivative  was  then  sub¬ 
stituted  into  (2)  and  solved,  yielding  a  simplified  radial 
profile  of  the  gradient  wind  (Fig.  9a).  While  modeling 
the  gradient  wind  as  a  quadratic  function  is  clearly  an 
oversimplification,  one  can  still  obtain  meaningful  esti¬ 
mates  of  super/subgradient  winds  in  Isabel’s  inner  core. 
From  Fig.  9a,  the  presence  of  supergradient  winds 
around  the  eyewall  region  ( r  >  20  km)  at  low  levels 
(—100  m)  is  consistent  with  the  inward  deceleration  of 
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Fig.  9.  (a)  Composite  wind  at  —100  m  (0-250-m  composite  level)  from  Hurricane  Isabel  on  13  Sep,  and  (b)  calculated  nonbalanced 
wind  in  a  cyclostrophic  vortex  with  specified  inflow  angle  such  that  the  radial  wind  converges  asymptotically  to  0.8  times  the  radius  of 
maximum  cyclostrophic  wind  (from  Willoughby  1990).  The  dotted  curve  represents  the  cyclostrophic,  or  gradient  wind;  the  solid  curve 
the  nonbalanced  tangential  wind;  the  shorter  dashed  curve  the  radial  wind;  and  the  longer  dashed  curve  the  difference  between  the 
balanced  and  nonbalanced  wind.  Wind  components  are  nondimensionalized  with  the  maximum  balanced  wind,  and  radius  is  nondi- 
mensionalized  with  the  radius  of  maximum  balanced  wind. 


the  radial  flow  (Smith  1980).  Although  the  simple  qua¬ 
dratic  function  falls  off  too  rapidly  at  outer  radii,  it 
suggests  that  the  boundary  layer  tangential  wind  was 
—15%  supergradient  near  the  RMW. 

The  transition  from  super-  to  subgradient  winds  oc¬ 
curs  at  — 20-km  radius  and  is  concurrent  with  a  reverse 
in  the  radial  acceleration  near  the  center  of  the  eye. 
This  pattern  is  similar  to  a  simple  hurricane  boundary 
layer  model  constructed  in  Willoughby  (1990),  shown 
in  Fig.  9b.  Willoughby’s  model  used  a  simple  log-spiral 
inflow  to  model  the  radial  wind,  but  the  predicted  wind 
structure  near  the  RMW  is  qualitatively  similar  to  the 
results  from  Isabel.  The  radial  acceleration  acts  effec¬ 
tively  like  an  extra  pressure,  decelerating  the  tangential 
wind  in  the  eye  and  accelerating  it  in  the  eyewall.  This 
analysis  is  also  qualitatively  similar  to  the  analytical 
model  presented  in  Kepert  (2001),  with  the  top  of  the 
inflow  layer  coinciding  with  the  location  of  maximum 
tangential  winds  and  a  supergradient  jet  below.  The 
results  are  similar  on  the  other  two  analysis  days  (not 
shown). 

b.  Thermal  wind  balance 

An  assessment  of  thermal  wind  balance  near  the 
RMW  at  the  boundary  layer  top  was  presented  in  M06 
for  13  September.  As  shown  in  section  4a,  the  assump¬ 
tions  made  in  that  paper  regarding  cyclostrophic  bal¬ 
ance  are  reasonably  satisfied  for  Isabel  above  the 
boundary  layer.  The  diagnostic  thermal  wind  equation 
relating  the  maximum  mean  tangential  wind  and  the 
radial  gradient  of  moist  entropy  at  the  top  of  the 


boundary  layer5  derived  by  Emanuel  [1986,  Eq.  (13); 
1997,  Eq.  (13)]  is  repeated  here  for  completeness: 

7  dSB 

V2m~~rm{TB-T0)  —  ,  (5) 

where  SB  is  the  moist  entropy  at  the  boundary  layer 
top,  /•  is  the  radius  from  the  vortex  center,  rm  is  the 
RMW,  Vm  is  the  maximum  mean  tangential  wind,  TB  is 
the  temperature  at  the  boundary  layer  top,  and  T0  is 
the  outflow  temperature. 

Similar  calculations  as  to  the  degree  of  thermal  wind 
balance  on  12  and  14  September  have  been  performed 
here.  The  results  from  all  three  days  are  summarized  in 
Table  2.  There  is  fair  agreement  between  the  observed 
quantities  on  the  13th  as  shown  in  M06.  On  the  14th, 
the  entropy  gradient  appears  too  weak  for  the  observed 
tangential  wind  and  outflow  temperature.  The  discrep¬ 
ancy  in  the  radial  0e  gradient  of  —1  K  (10  km)-1  is 
within  the  range  of  uncertainty  for  this  measurement 
and  the  sensitivity  to  this  slight  discrepancy  on  this  day 
yields  either  a  30-K  decrease  in  outflow  temperature,  or 
a  10  m  s-1  decrease  in  the  tangential  wind  for  exact 
axisymmetric  thermal  wind  balance.  A  more  distinct 


5  This  theoretical  thermal  wind  relationship  is  strictly  valid  at 
the  height  of  the  maximum  wind  at  the  RMW.  From  Fig.  4.  the  top 
of  the  inflow  layer  and  the  height  of  the  maximum  wind  are  at 
approximately  the  same  level  near  the  RMW  on  all  three  days. 
This  altitude  is  therefore  defined  as  the  boundary  layer  top  fol¬ 
lowing  Montgomery  et  al.  (2001).  Mixed-layer  depths  for  virtual 
potential  temperature  were  much  shallower  (—150  m),  however, 
as  shown  in  M06. 
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Table  2.  Results  of  thermal  wind  diagnostic  from  Eq.  (5)  for  12-14  Sep.  Top  row  and  first  column  values  are  observed  quantities. 
Cell  contents  in  boldface  show  the  results  of  calculating  for  the  remaining  parameter  assuming  the  other  two  quantities  are  known. 


12  Sep,  25-km  RMW 
Vm=  80  m  s"1 

Ta  =  -65°C 

13  Sep,  45 -km  RMW 
V,„=  76  m  s"1 

T0  =  -58°C 


Ta  =  -65°C 

d0e/dr  =  -10.2  K  (10  km)"1 


Ta  =  -58°C 

d0e/dr  =  -5.7  K  (10  km)"1 


dd^dr  =  -5  K  (10  km)  1 

Ta  =  —158°C 
Vm  =  56  ms"1 

dd^dr  =  -6  K  (10  km)"1 

Ta  =  —  57°C 
Vm  =  74  ms"1 


14  Sep,  50-km  RMW 
Vm=  74  m  s"1 
Ta  =  -56°C 


T0  =  -56°C 

dOe/dr  =  -4.7  K  (10  km)"1 


dd/dr  =  -3.5  K  (10  km)"1 

T0  =  —  83°C 
Vm  =  64  ms"1 


discrepancy  between  the  wind  and  entropy  fields  at  the 
RMW  on  12  September  is  evident.  Because  of  uncer¬ 
tainties  in  the  inferred  radial  entropy  gradient  and  cal¬ 
culated  outflow  temperature  (described  further  in  sec¬ 
tion  5b),  it  is  difficult  to  assert  with  confidence  that  the 
storm  was  significantly  unbalanced  in  this  regard.  The 
steeper  radial  entropy  gradient  near  25-km  radius,  how¬ 
ever,  and  the  apparent  consistency  at  the  larger  RMW 
(45  km)  on  the  13th  suggests  that  some  adjustment  to¬ 
ward  thermal  wind  balance  may  have  taken  place. 
However,  with  the  limited  temporal  continuity,  the  use 
of  pseudoequivalent  potential  temperature,  and  the 
above-mentioned  uncertainties,  further  speculation  on 
the  matter  is  not  justified. 

5.  Potential  intensity  analysis 

A  detailed  PI  analysis  of  Hurricane  Isabel  on  13  Sep¬ 
tember  was  performed  in  M06,  where  it  was  argued  that 
Isabel  exceeded  the  EPI  a  priori  maximum  mean  tan¬ 
gential  wind  within  the  range  of  10-35  m  s-1  on  13 
September.  The  range  of  values  depends  on  the  degree 


of  negative  ocean  feedback,  the  numerical  value  of  en¬ 
thalpy  and  momentum  exchange  coefficients  at  high 
surface  wind  speeds,  and  dissipative  heating.  Here  we 
continue  that  analysis  with  available  data  from  12  and 
14  September.  For  an  EPI  prediction  of  maximum 
mean  tangential  wind  at  the  top  of  the  boundary  layer, 
one  needs  estimates  of  1)  the  SST  underneath  the  hur¬ 
ricane  core,  2)  the  upper-level  exhaust  (outflow)  tem¬ 
perature  where  parcels  undergo  radiational  cooling  to 
space,  3)  the  ambient  RH  near  the  sea  surface,  and  4) 
the  ratio  of  bulk  enthalpy  to  momentum  exchange  co¬ 
efficients  at  high  wind  speeds.  A  constant  RH  of  80% 
was  used  in  this  study  as  a  simple  expedient  [see  Table 
3  for  ambient  RH  measurements,  and  M06  and  Camp 
and  Montgomery  (2001)  for  a  discussion  of  this  assump¬ 
tion].  The  other  parameters  are  examined  in  the  fol¬ 
lowing  sections. 

a.  Sea  surface  temperature 

The  average  SST  over  the  week  prior  to  Isabel’s  pas¬ 
sage  derived  from  Advanced  Very  High  Resolution  Ra- 


Table  3.  Observed  and  calculated  environmental  temperatures  and  RH  for  Hurricane  Isabel  from  12  to  14  Sep  2003.  RH  is  the 
average  in  the  lowest  500  m.  Temperature  at  the  boundary  layer  top  ( TB ),  radial-wind-weighted  outflow  temperature  (Ta  RW), 
equilibrium-level  outflow  temperature  ( Ta  EL),  f/ln(0c)-weighted  outflow  temperature  ( Ta  EPI),  and  de  at  the  outflow  layer  top  and 
bottom  are  also  shown. 


Radial  outflow  layer 


from  Altitude 


Sounding 

center 

(km) 

RH 

(%) 

tb 

(°C) 

Ta  RW 
(°C) 

Ta  EL 
(°C) 

Ta  EPI 
(°C) 

0e  Bottom 
(K) 

bottom 

(m) 

6e  top 
(K) 

Altitude 
top  (m) 

TJSJ  2400:00  UTC  12  Sep  (SW) 

861 

82 

21.4 

-54.5 

-58 

-64.6 

351.7 

12  420 

363.0 

15  030 

G-IV  2227:14  UTC  13  Sep  (NE) 

341 

72.6 

24.9 

-54.2 

-59 

-52.9 

349.7 

11  140 

358.8 

14  060 

G-IV  2004:13  UTC  13  Sep  (NW) 

525 

82.5 

24.5 

-43.4 

-44 

-54.4 

349.7 

11700 

354.9 

13  210 

G-IV  0033:36  UTC  13  Sep  (SE) 

347 

77.7 

23.8 

-55.7 

-50 

-56.5 

349.8 

12  090 

359.2 

14  359 

TJSJ  2400:00  UTC  13  Sep  (SW) 

584 

84.2 

24.3 

-70.7 

-78 

-79.5 

349.7 

14  109 

359.2 

16  539 

Average  of  7  soundings  on  13  Sep 

653.8 

80.7 

24.1 

-52.9 

-56.7 

-57.8 

349.7 

12  389 

358.0 

14  801 

G-IV  2004:42  UTC  14  Sep  (N) 

607 

81.0 

24.3 

-51.2 

-47 

-56.2 

349.8 

11  860 

357.8 

13  470 

G-IV  1921:57  UTC  14  Sep  (NE) 

371 

83.9 

24.0 

-45.8 

-52 

-51.3 

349.7 

11  520 

357.5 

13  290 

G-IV  1836:44  UTC  14  Sep  (SE) 

726 

72.9 

24.4 

-47.1 

-48 

-61.5 

349.7 

12  280 

359.3 

13  770 

TJSJ  2400:00  14  Sep  (S) 

726 

71.8 

23.2 

-39 

-61.2 

-49.9 

349.7 

11  079 

359.3 

14  009 

Average  of  7  soundings  on  14  Sep 

609.7 

81.9 

24 

-50.3 

-49.4 

-55.5 

349.8 

11811 

356.8 

13  417 
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Sea  Surface  Temperature  (°C) 


Fig.  10.  SST  derived  from  AVHHR  satellite  (average  SST  over  4-10  Sep  in  color)  and  NOAA  WP-3D  downward-pointing  radi¬ 
ometer  (thin  line,  from  —1800  UTC  13  and  14  Sep).  Tracks  of  Hurricanes  Fabian  (dashed  best  track,  from  2  to  5  Sep)  and  Isabel  (dashed 
best  track,  with  thick  white,  solid  line  indicating  analysis  periods  from  1600  to  2300  UTC  12-14  Sep)  are  shown  for  reference. 


diometer  (AVHRR)  satellite  measurements  (McMillin 
and  Crosby  1984)  is  shown  in  Fig.  10  with  the  analysis 
periods  highlighted  in  white.  The  prestorm  satellite- 
derived  temperatures  are  well  validated  by  airborne  ra¬ 
diometer  measurements  (Moss  1978)  recorded  on 
board  the  P-3  NOAA  43  aircraft  at  —60  and  —120  m 
above  sea  level  to  the  northwest  (in  advance)  of  the 
storm  center  on  14  September.  Radiometer  measure¬ 
ments  taken  in  the  rear  of  Isabel  on  13  September  mea¬ 
sured  a  reduction  in  SST  by  —  1°-2°C  relative  to  satel¬ 
lite-derived  SST  estimates  prior  to  Isabel's  passage.  As 
discussed  in  M06,  this  reduction  is  believed  to  be 
caused  by  shear-induced  turbulent  upwelling  (Price 
1981;  Shay  and  Elsberry  1987;  Emanuel  et  al.  2004),  an 
effect  that  is  not  included  in  current  EPI  theory.  Isabel 
encountered  the  cool  wake  of  Fabian  in  the  central 
region  of  the  plot  during  the  IOP  on  the  13th,  and  a 
constant  SST  of  27.5°C  was  employed  in  M06  to  repre¬ 
sent  the  average  SST  for  the  a  priori  PI  estimates.  SSTs 
on  the  12th  had  values  near  the  inner  core  of  — 28.5°C, 
and  Isabel  began  to  encounter  slightly  warmer  water  on 
the  14th,  with  the  estimated  SST  near  — 29°C. 

M06  presented  arguments  that  the  maximum  sus¬ 
tained  tangential  wind  rapidly  adjusts  to  SST  changes 
underneath  the  eyewall,  consistent  with  EPI  theory.  A 
numerical  modeling  experiment  was  performed  to  test 
the  effect  of  a  cool  wake,  similar  to  that  experienced  by 
Isabel.  This  experiment,  summarized  in  the  appendix, 


suggests  that  the  intensity  response  is  on  the  order  of  a 
few  hours,  and  it  is  therefore  likely  that  Isabel’s  Fmax 
rapidly  adjusted  to  the  lower  SST. 

b.  Outflow  temperature 

The  outflow  temperature  was  calculated  in  M06  in 
three  different  ways:  as  a  radial-wind-weighted  tem¬ 
perature  across  the  storm  outflow  (PM03);  as  an  equi¬ 
librium-level  temperature  (the  temperature  at  which  a 
virtual  parcel  starting  from  ambient  surface  state 
achieves  the  same  environmental  temperature  after  lift¬ 
ing  by  pseudoadiabatic  ascent);  and  as  a  d\n(0e)- 
weighted  temperature  following  the  original  theoretical 
definition  (Emanuel  1986;  see  M06  for  more  details 
relevant  to  this  study).  The  first  approach  provides  per¬ 
haps  the  most  empirical  of  the  three  calculations,  in  that 
it  is  a  direct  measurement  of  the  temperature  in  the 
radial  outflow  of  the  storm.  Its  deficiency  with  obser¬ 
vational  data  is  that  the  outflow  is  asymmetric  and  of¬ 
ten  concentrated  in  jets  (see  Fig.  6  of  M06;  Vladimirov 
et  al.  2001);  thus  uneven  sampling  may  bias  the  esti¬ 
mate.  The  second  provides  a  plausible  proxy,  in  that 
parcels  that  rapidly  ascend  in  the  eyewall  must  return  to 
a  level  of  neutral  buoyancy  as  they  recede  from  the 
storm,  cooling  gradually  and  sinking  at  large  distances. 
The  equilibrium  level  of  a  particular  sounding  depends 
heavily  on  the  surface  characteristics;  however,  it  may 
be  biased  if  these  characteristics  are  significantly  differ- 
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(a)  GOES  IR  Image  and  Soundings,  12  September 


(b)  14  September 


T800Z  13 L  ISABEL 


W7T5/03  00002  13L  ISABEL 


126Z  TRMM  over; 
045Z  GOES-12  il 


Naval  Research  Laboratory 


navy.mil/sat_products.htinl  Naval  Research  Laboratoi 


Fig.  11.  GOES  infrared  satellite  imagery  (courtesy  of  NRL/Monterey)  at  (a)  2045  UTC  12  Sep  and  (b)  2145  UTC  14  Sep  showing 
extent  of  hurricane  outflow,  and  the  dropwindsonde  profile  locations  (red  dots  with  UTC  time  labels)  used  in  the  outflow  temperature 
calculation. 


ent  than  that  found  in  the  hurricane  environment  (i.e., 
if  the  sounding  was  taken  over  land  instead  of  ocean). 
The  third  definition  emerges  naturally  in  the  EPI  for¬ 
mulation,  and  is  based  on  the  “cool  reservoir”  of  the 
Carnot  engine  concept.  It  is  the  most  relevant  calcula¬ 
tion  and  forms  the  primary  basis  for  the  outflow  tem¬ 
peratures  presented  here. 

The  average  of  seven  independently  calculated  out¬ 
flow  temperatures  on  13  September  was  —  58°C  and 
was  used  as  the  primary  system-scale  outflow  tempera¬ 
ture  for  that  day  (see  Fig.  6  of  M06).  Since  the  NOAA 
G-IV  jet  was  not  deployed  on  12  September,  a  single 
radiosonde  (TJSJ)  from  San  Juan,  Puerto  Rico,  pro¬ 
vides  the  only  in  situ  estimate  of  the  outflow  tempera¬ 
ture  on  that  day.  The  San  Juan  sounding  appeared  to  be 
an  outlier  on  the  13th,  possibly  due  to  the  high  tropo- 
pause  height  at  low  latitude,  the  high  equilibrium  level 
since  the  rawinsonde  was  launched  from  land,  and/or 
the  different  sensor  type  from  the  GPS  dropwind- 
sondes.  The  calculated  outflow  temperatures  from  this 
sounding  are  shown  in  the  first  row  of  Table  3.  In  this 
case,  the  radial  wind  weighted  outflow  (RW)  and  equi¬ 
librium  level  (EL)  temperatures  are  warmer  than  the 
dln(9e)  weighted  temperature  (EPI).  The  EPI  tempera¬ 
ture  of  — 65°C  is  reasonable,  however,  yielding  a  colder 
outflow  than  on  the  13th  or  14th  due  to  higher  6e  inte¬ 
gral  limits  (columns  7  and  9  in  Table  3)  determined  by 
the  composite  mean  6e  in  the  eyewall.  Infrared  satellite 


imagery  provides  a  qualitative  validation  of  this  esti¬ 
mate,  showing  generally  colder  cloud-top  temperatures 
than  on  the  13th  or  14th,  but  no  pixels  below  — 75°C 
(Fig.  11).  This  single  estimate  results  in  an  additional 
uncertainty  in  the  a  priori  maximum  mean  tangential 
wind  of  3-6  m  s_1,  depending  on  the  specific  assump¬ 
tions  used  in  the  calculation. 

On  14  September,  G-IV  reconnaissance  provides  a 
comprehensive  look  at  the  ambient  environment,  al¬ 
lowing  for  a  more  robust  calculation  of  the  outflow 
temperature  similar  to  13  September.  The  bottom  rows 
of  Table  3  show  a  subset  of  the  seven  soundings  used  to 
calculate  the  outflow  on  this  day.  The  three  methods 
described  above  yield  similar  estimates  near  — 56°C. 
This  also  consistent  with  the  value  calculated  on  the 
13th  when  the  storm  had  similar  thermodynamic  struc¬ 
ture  at  the  eyewall,  and  with  the  apparent  decrease  in 
cloud-top  heights  as  seen  in  Fig.  11. 

c.  Bulk  enthalpy  and  momentum  exchange 
coefficients 

The  complex  nature  and  behavior  of  the  ocean  sur¬ 
face  at  high  wind  speeds  is  one  of  the  least  understood 
aspects  of  hurricane  science.  Observations  of  turbulent 
fluxes  of  enthalpy  and  momentum  in  the  hurricane 
boundary  layer  at  the  highest  wind  speeds  to  date  (—30 
m  s-1)  were  collected  during  the  2002-04  CBLAST 
field  program.  These  fluxes  were  then  used  to  calculate 
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Fig.  12.  Photograph  of  sea  surface  in  Hurricane  Isabel  taken  from  NOAA  43  during  stepped  descent  pattern  below  400-m  altitude. 

Surface  wind  speed  is  20-30  m  s-1.  (Photo  courtesy  of  M.  Black.) 


bulk  surface  exchange  coefficients  as  a  way  of  param¬ 
eterizing  the  air-sea  interaction.  These  results  suggest 
that  current  estimates  of  the  ratio  of  exchange  coeffi¬ 
cients  for  minimal  hurricane-force  wind  speeds  are  ap¬ 
proximately  0.7  (Black  et  al.  2007;  Drennan  et  al.  2007; 
French  et  al.  2007).  Direct  measurements  were  not  ob¬ 
tained  above  these  wind  speeds,  however,  because  of 
the  inherent  difficulty  and  danger  in  collecting  airborne 
near-surface  measurements  at  higher  wind  speeds.  A 
photograph  of  the  sea  surface  at  hurricane-force  winds 
taken  during  the  low-level  (~60-m  altitude)  passes  on 
13  September  is  shown  in  Fig.  12.  Wind  streaks  and 
waves  of  varying  size  and  age  make  parameterizing  this 
complex  air-sea  interaction  challenging,  to  say  the 
least. 

At  even  higher  wind  speeds,  dropsonde  observations 
(Powell  et  al.  2003)  and  laboratory  tank  experiments  by 
Donelan  et  al.  (2004)  suggest  that  the  aerodynamic 
roughness  of  the  surface,  and  therefore  CD,  reaches  a 
limiting  value.  Direct  observations  of  the  enthalpy  flux 


at  extreme  wind  speeds  have  not  yet  been  made,  but 
alternative  means  to  determine  this  coefficient  via  a 
budget  residual  method  have  been  performed  using 
flight-level  data  (e.g.,  Hawkins  and  Imbembo  1976)  and 
are  currently  underway  including  dropwindsonde  data 
(K.  Emanuel  2005,  personal  communication).  The  ef¬ 
fects  of  sea  spray  on  surface  fluxes  at  wind  speeds 
above  30  m  s-1  are  still  poorly  understood,  but  likely 
play  an  important  role  in  the  air-sea  interaction  in  ma¬ 
jor  hurricanes  (Fairall  et  al.  2003).  Given  these  uncer¬ 
tainties,  the  ratio  of  bulk  enthalpy  and  momentum  ex¬ 
change  coefficients  is  assumed  to  be  unity  ( CKICD  =  1) 
for  the  both  the  upper  bound  and  primary  PI  estimates 
in  this  paper,  with  a  lower  bound  of  0.5;  this  range 
yields  an  uncertainty  in  PI  of  —20  m  s-1. 

d.  PI  estimates:  Azimuthal  mean  Vmax  at  boundary 
layer  top 

Figure  13  shows  the  predicted  Fmax  from  the  a  priori 
EPI  theory  for  varying  outflow  temperatures  and  near- 
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core  SSTs  with  a  constant  RH  =  80%.  All  calculations 
were  performed  using  the  same  assumptions  as  in  M06, 
where  the  upper  bound  for  the  predicted  Fmax  uses  a 
CKICD  value  of  one,  assumes  a  full  contribution  from 
dissipative  heating  (Bister  and  Emanuel  1998),  and  as¬ 
sumes  no  storm-induced  ocean  cooling  (Emanuel  et  al. 
2004).  The  lower  bound  assumes  a  CKICD  equal  to  0.5, 
and  that  dissipative  heating  is  entirely  offset  by  the 
ocean  cooling  feedback.  The  primary  estimate  compro¬ 
mises  these  assumptions,  with  a  bulk  exchange  coeffi¬ 
cient  ratio  of  one,  and  offsetting  dissipative  heating  and 
ocean  cooling.  This  yields  a  theoretically  predicted  56.6 
m  s-1  EPI  for  the  environmental  conditions  on  13  Sep¬ 
tember  (Fig.  13b),  resulting  in  a  10-35  m  s  1  discrep¬ 
ancy  between  theory  and  observations  as  Isabel  crossed 
the  cool  wake  of  Fabian. 

Using  these  same  assumptions  one  obtains  a  primary 
theoretical  PI  estimate  of  61.2  m  s-1  for  12  September. 
This  value  is  —19  m  s"1  below  the  observed  mean  tan¬ 
gential  wind  at  the  top  of  the  boundary  layer  of  80 
m  s-1.  The  upper  and  lower  bound  estimates  yield  a 
range  of  42-76  m  s_1,  including  the  outflow  tempera¬ 
ture  uncertainty.  The  highest  value  is  within  the  stan¬ 
dard  deviation  (10  m  s-1)  of  the  mean  tangential  wind 
estimate  in  this  case.  As  was  shown  in  section  4,  the 
boundary  layer  9e  structure  on  12  September  conforms 
more  closely  to  a  more  classic  (in-up-out)  thermody¬ 
namic  pathway,  where  substantial  6e  gain  is  achieved  in 
the  outer  core  and  underneath  the  eyewall.  Additional 
high  de  in  the  eye  and  the  presence  of  eye  and  eyewall 
mesovortices  (Fig.  3)  may  be  partially  responsible  for 
the  storm  exceeding  the  a  priori  PI,  however. 

Warmer  sea  surface  temperatures  yield  slightly 
higher  PI  values  on  14  September  than  those  on  13 
September,  but  weaker  than  that  found  on  12  Septem¬ 
ber.  Using  the  same  assumptions  for  the  upper  and 
lower  bounds  yields  a  range  of  42-69  m  s-1,  with  a  best 
estimate  of  59.0  m  s-1.  The  PI  is  therefore  exceeded  on 
this  day  as  well,  with  the  maximum  value  of  Qe  in  the 
eye  diminished,  but  the  overall  area  of  the  potential 
entropy  reservoir  was  large,  similar  to  the  previous  day. 
While  the  size  of  the  eye  likely  plays  a  role  in  the  en¬ 
ergetics  of  the  system,  it  is  unknown  at  this  point  how 
much  of  that  energy  can  be  tapped  by  the  hurricane. 


<— 


Fio.  13.  Theoretically  predicted  azimuthal  mean  Vmax  at  the 
boundary  layer  top  for  varying  outflow  temperature  and  near- 
core  SST  with  a  constant  RH  =  80%.  The  “X”  indicates  the 
primary  potential  intensity  estimate  for  the  observed  near  envi¬ 
ronment  around  Isabel.  The  dark  solid  curve  represents  the  av¬ 


erage  storm-relative  tangential  wind  speed  at  the  top  of  the 
boundary  layer  derived  from  the  dropwindsonde  measurements. 
The  shading  represents  the  standard  deviation  of  this  mean  value. 
The  a  priori  EPI  estimates  assume  CKICD  =  1.  See  text  for  as¬ 
sumptions  leading  to  upper  and  lower  bound  estimates. 
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6.  Conclusions 

An  analysis  of  the  structure,  evolution,  and  PI  of 
Hurricane  Isabel  near  its  maximum  intensity  was  per¬ 
formed  using  a  combination  of  dropwindsonde,  flight- 
level,  satellite,  and  airborne  radar  data.  The  analysis 
suggests  a  nearly  steady-state  intensity,  with  an  expand¬ 
ing  tangential  wind  field  accompanied  by  dramatic 
changes  in  the  secondary  circulation  and  moist  entropy 
structure.  A  distinct  increase  in  near-surface  radial  in¬ 
flow,  eyewall  vertical  velocity,  and  low-level  0e  in  the 
eye  from  12  to  13  September  suggests  a  buildup  of 
moist  entropy  due  to  latent  heat  fluxes  in  the  relatively 
quiescent,  low-pressure  eye,  which  is  then  accessed  by 
parcels  that  are  able  to  penetrate  the  eyewall,  via  the 
intense  low-level  mean  inflow  or  asymmetric  mixing 
associated  with  mesovortices.  The  injection  of  this  heat 
energy  into  the  eyewall  supports  the  strengthening  of 
the  secondary  circulation,  both  in  an  axisymmetric 
sense  (Shapiro  and  Willoughby  1982;  PM03)  and  as  a 
local  buoyancy  source  (Braun  2002;  Eastin  et  al.  2005), 
resulting  in,  or  possibly  concurrent  with,  increased  ad- 
vection  of  higher  momentum  air  from  the  outer  core 
and  expansion  of  the  wind  field.  Despite  a  cooler  SST 
and  a  rise  in  central  pressure,  Isabel  maintained  76 
ms  1  mean  tangential  winds  on  13  September.  The 
storm  weakened  slightly  on  14  September,  and  exhib¬ 
ited  a  flattening  of  the  radial  entropy  gradient  in  the 
eyewall  region  consistent  with  mixing  across  the  eye- 
eyewall  interface. 

The  inner-core  Rossby-Ertel  potential  vorticity 
structure  obtained  through  the  Barnes  analysis  pro¬ 
vides  a  new  look  at  the  details  of  this  important  dy¬ 
namical  quantity  in  a  category  5  hurricane.  Significant 
absolute  vertical  vorticity  (>6  X  10-3  s-1)  and  PV  (>60 
PVU)  are  found  radially  inward  of  the  eyewall  on  each 
day,  with  a  pronounced  ring  structure  typically  associ¬ 
ated  with  very  intense  storms.  This  ring,  or  bowl,  of 
high  PV  supports  the  necessary  condition  for  combined 
barotropic/baroclinic  instability  (Montgomery  and  Sha¬ 
piro  1995),  which  has  been  shown  to  be  a  probable 
mechanism  for  producing  hurricane  mesovortices  and 
polygonal  eyewalls  (Schubert  et  al.  1999;  Montgomery 
et  al.  2002;  Kossin  and  Schubert  2004).  This  breakdown 
of  the  vortex  sheet  allows  for  mixing  at  the  eye-eyewall 
interface  in  conjunction  with  the  frictionally  driven 
mean  inflow,  providing  a  plausible  dynamical  pathway 
for  the  thermodynamic  augmentation  of  the  hurricane 
engine  described  in  this  paper.  A  secondary  PV  maxi¬ 
mum  is  found  near  3-km  altitude  in  the  eye,  coupled 
with  a  thermal  inversion  and  decrease  in  relative  hu¬ 
midity.  This  feature  has  been  identified  in  recent  nu¬ 
merical  simulations  of  intense  hurricanes  (Chen  and 


Yau  2001;  Wang  and  Zhang  2003),  but  its  dynamical 
significance  has  not  been  explored  thoroughly  to  the 
authors’  knowledge.  This  remains  a  topic  for  further 
research. 

An  analysis  of  gradient  wind  balance  on  13  Septem¬ 
ber  suggests  that  the  winds  in  the  boundary  layer  were 
supergradient  in  the  eyewall  region,  but  transitioned  to 
subgradient  inside  the  eye.  This  is  consistent  with  the 
radial  accelerations  observed  in  the  axisymmetric  wind 
composite,  and  the  inferred  radial  structure  is  similar  to 
the  simplified  log-spiral  inflow  boundary  layer  model  of 
Willoughby  (1990).  At  ~100-m  altitude,  the  contribu¬ 
tion  from  the  radial  advection  and  eddy/friction  terms 
was  significant,  opposing  the  cyclostrophic  radial  pres¬ 
sure  gradient  and  reducing  the  central  pressure  deficit 
by  —15  hPa.  At  2-km  altitude,  the  observed  radial  pres¬ 
sure  gradient  and  integrated  pressure  deficit  agree  well, 
indicating  that  the  storm  was  near  gradient  balance 
above  the  boundary  layer. 

Thermal  wind  balance  was  generally  difficult  to  as¬ 
sess  using  solely  the  dropwindsonde  data  because  of  the 
diagnostic  equation’s  sensitivity  to  the  radial  gradient 
of  moist  entropy  in  the  eyewall.  Results  suggested  that 
the  observed  moist  entropy  gradient  and  outflow  tem¬ 
perature  were  roughly  consistent  with  the  maximum 
mean  tangential  wind  on  13  September,  but  less  so  on 
14  September.  The  analysis  also  revealed  a  distinct  dis¬ 
crepancy  on  12  September  in  the  observed  radial  en¬ 
tropy  gradient  from  that  calculated  by  the  thermal  wind 
relationship.  The  steeper  radial  entropy  gradient  at  the 
former  RMW  (25  km)  and  the  apparent  balance  at  the 
larger  RMW  (45  km)  by  the  next  day  suggest  that  some 
adjustment  toward  thermal  wind  balance  may  have 
taken  place,  but  limited  temporal  continuity  and  uncer¬ 
tainty  in  the  entropy  gradient  and  outflow  temperature 
preclude  further  speculation  as  to  the  details  of  this 
adjustment  process. 

We  have  focused  this  paper  on  comparisons  with  EPI 
theory  since  air-sea  heat  and  momentum  fluxes  are  not 
explicitly  included  in  Miller’s  (1958)  or  Holland’s 
(1997)  PI  theories,  nor  do  these  theories  offer  predic¬ 
tions  for  Vmax.  This  study  suggests  that  Hurricane  Isa¬ 
bel’s  structure  and  intensity  was  largely  consistent  with 
many  of  the  central  concepts  of  the  EPI  theory  regard¬ 
ing  boundary  layer  balance,  and  the  hurricane  as  a  Car¬ 
not  engine.  However,  the  evidence  that  Isabel  exceeded 
Emanuel’s  PI  formulation  indicates  that  revisions  to 
this  important  dynamical  concept  need  to  be  addressed. 
This  study  provides  further  evidence  that  entropy  pro¬ 
duction  inside  the  RMW  in  the  low-level  eye  provides 
an  active  role  in  hurricane  intensification,  serving  as  an 
important  and  overlooked  energy  source  for  the  hurri¬ 
cane  engine,  and  challenges  classic  theories  about  the 
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degree  of  moist  entropy  obtained  from  the  ocean  in  the 
outer  core  and  underneath  the  eyewall. 

While  the  conclusions  presented  here  appear  robust 
for  this  dataset,  further  research  with  additional  case 
studies  and  numerical  modeling  is  needed  to  elucidate 
the  various  thermodynamic  pathways  available  to  the 
tropical  cyclone  for  maintenance  and  growth.  The  rela¬ 
tive  importance  of  the  size  of  the  eye,  the  amount  of  eye 
6e  that  is  utilized  by  the  hurricane  via  injection  into  the 
eyewall  cloud,  and  the  how  important  the  mechanisms 
described  in  this  paper  may  be  in  weaker  storms  are 
interesting  and  open  questions  for  future  research. 
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APPENDIX 

The  Fmax  Response  to  SST  Change 

The  scaling  argument  offered  in  footnote  1  of  M06 
and  dynamical  considerations  based  on  the  intrinsic  at¬ 
mospheric  time  scales  near  the  RMW  of  an  intense 
vortex  would  both  suggest  that  the  adjustment  time  for 
Fmax  due  to  a  reduced  SST  around  the  eyewall  region 
should  be  rapid  (on  the  order  of  a  few  hours).  To  test 
this  hypothesis  we  asked  Dr.  John  Persing  to  assist  us 
by  performing  two  experiments  using  the  axisymmetric 
Rotunno  and  Emanuel  (1987)  numerical  model  at  the 
"4x”  resolution  detailed  in  PM03  with  a  modified 
sounding.  The  two  experiments  were  a  control  run 
(CTL)  using  an  SST  of  26.1°C  and  a  sensitivity  experi¬ 
ment  (EXP)  with  an  SST  lowered  by  1°C  within  60  km 
from  the  vortex  center  after  6.5  simulation  days  when 
the  control  experiment  had  attained  a  quasi-steady  tan¬ 
gential  velocity  time  series  (not  shown).  The  120-km- 
diameter  disk  of  lower  SST  is  a  simple  way  of  mimick- 
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Fig.  Al.  Difference  in  the  maximum  tangential  wind  (Vmax) 
between  the  experimental  run  (EXP)  and  control  run  (CTL)  over 
time  of  the  axisymmetric  hurricane  model  described  in  PM03  and 
Rotunno  and  Emanuel  (1987).  The  EXP  run  imposes  1°C  SST 
reduction  underneath  the  eyewall  at  6.5  days  into  the  simulation 
to  mimic  the  passage  over  a  cool  hurricane  wake. 

ing  the  passage  of  Isabel  over  Fabian’s  wake  during  the 
observation  period. 

The  outcome  of  these  two  experiments  is  summa¬ 
rized  in  Fig.  Al,  which  displays  the  difference  between 
the  control  and  sensitivity  experiment  of  the  maximum 
tangential  velocity  (in  meters  per  second)  as  a  function 
of  time  (in  days).  As  expected,  before  6.5  days  simula¬ 
tion  time  the  difference  is  identically  zero.  Immediately 
after  the  cooling  ring  is  imposed  in  the  model,  however, 
the  difference  plot  shows  a  rapid  and  nearly  linear  de¬ 
crease  in  Fmax  of  several  meters  per  second  in  approxi¬ 
mately  2  h.  EPI  theory  predicts  an  approximately  2 
ms^1  decrease  per  1°C  decrease  near  the  RMW  (see 
Fig.  Al).  Although  this  simulated  temporary  weaken¬ 
ing  is  stronger  than  predicted  using  steady-state  theory, 
the  main  point  here  is  that  Fmax  adjusts  nearly  instan¬ 
taneously  in  the  direction  of  the  imposed  SST  change. 

On  longer  time  scales,  the  intrinsic  variability  be¬ 
tween  the  two  runs  does  manifest  itself  and  it  is  not  yet 
clear  to  us  how  much  of  this  is  due  to  the  intrinsic  chaos 
of  an  intense  hurricane  vortex  and  model  numerics. 
Although  the  difference  plot  does  show  that  the  vortex 
temporarily  rebounds  from  this  local  decrease  in  SST 
on  a  time  scale  of  1.5  days,  only  the  early  time  tendency 
is  pertinent  to  the  crossing  of  Fabian’s  wake.  These 
numerical  simulations  support  the  hypothesis  that  Fmax 
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responds  to  the  lowered  SST  surrounding  the  eyewall 

region  on  a  relatively  short  (hourly)  time  scale  and  that 

Vmax  is  not  controlled  by  the  SST  in  the  outer  region  of 

the  storm. 
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